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'       Abstract 

Variational  properties  for  the  spectral  radius  and  spectrjd  abscissa  of  an  analytic 
matrix  valued  mapping  A:C  >-♦  €"*"  are  considered.  A  notion  of  directional  differen- 
tiability is  introduced  that  allows  us  to  exploit  the  perturbation  results  of  Part  1  of  this 
study.  Lower  bounds  for  the  directional  derivative  are  established  which  yield  formu- 
las for  the  directional  derivative  when  a  natural  nondegeneracy  condition  is  satisfied. 
These  formulas  are  interpreted  in  the  extreme  cases  where  the  eigenvalues  attaining 
either  the  spectral  radius  or  the  spectral  abscissa  are  nonderogatory  or  semisimple 
(nondefective).  We  conclude  by  investigating  the  relationship  with  the  proximal  nor- 
mal subdifferential.  In  particular,  it  is  shown  that  the  proximal  normjJ  subdifferential 
is  always  nonempty. 

1      Introduction 

In  this  article,  we  develop  a  few  of  the  variational  properties  of  the  spectral  radius,  p,  and  the 
spectral  abscissa,  a,  for  the  analytic  matrix  valued  mapping  A:C'  >-*  C"^".  The  mappings 
p  and  Q  are  given  by  the  formulas 

p(;):=max{|A|:AeS(r)}, 

and 

q{z)  :=max{ReA  :  A  €  E(z)}, 

where  the  multifunction  E:  (f  t-^  (C">^"  is  the  spectrum  of  A(z),  i.e.  the  roots  of  the 
characteristic  polynomial 

P(z,X)  =  det(XI-A(z)). 
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As  discussed  in  Part  I  [3],  the  mappings  a  and  p,  although  continuous,  are  not,  in  general, 
lipschitzian  near  values  of  z  for  which  P  has  multiple  roots.  The  nonlipschitzian  behav- 
ior arises  from  the  bifurcation  of  these  multiple  roots.  Along  analytic  curves  in  C  these 
roots  are  algebroidal  functions  and  as  such  can  be  described  as  the  branches  of  an  analytic 
function  with  each  branch  being  representable  as  a  power  series  in  fractional  powers  of  the 
perturbation.  It  is  these  classical  results  that  inspire  us  to  introduce  a  notion  of  directional 
differentiability  which  depends  in  an  essential  way  on  the  underlying  analyticity  of  ^. 
For  /:(["  H-IR,  define  f''(z\  ■):€'>-*  MU  {±00}  by 

f'iz;d):=      mf     lim  mf  ^^^M^lM,  (j) 

where 

r(z,d):=  {r.€<-^C'    :   7    is  analytic  ,  7(0)  =  r,  and  7'(0)  =  d}.  (2) 

The  notation  €  |  0  implies  that  only  real  (and  positive)  values  of  t  are  to  be  considered. 
The  superscript  h  in  (1)  is  used  to  emphasize  that  f'^(z\d)  depends  only  on  the  analytic 
(holomorphic)  curves  in  C  that  pass  through  z.  At  first  glance  it  may  appear  that  this 
notion  of  directional  differentiability  is  not  sufficiently  robust  to  capture  the  full  range  of 
first  order  variational  behavior  that  is  possible  for  nonlipschitzian  functions.  However,  as 
we  hope  to  demonstrate,  it  is  sufficiently  robust  for  the  analysis  of  the  functions  p  and  q.  In 
order  to  see  why  this  is  so,  we  need  to  recall  certain  properties  of  the  Puiseux-Newton  series 
expansions  for  the  roots  of  P,  which  were  introduced  by  way  of  example  in  the  discussion 
following  Lemma  1  in  Part  I  of  this  study.  Here  we  provide  a  little  more  detail. 

Let  z,  d  be  in  C  and  choose  7  G  r(z,d).   The  roots  of  P(7(f),A)  are  given  by  series 
expansions  in  fractional  powers  of  e,  that  is,  Puiseux-Newton  series  of  the  form 

A(f)  :=Ao +  «''  +  •■■,  (3) 

where  p  :=  q/r  >  1/n  for  some  pair  of  relatively  prime  integers  q  and  r.  Associated  with 
each  such  value  of  p  there  are  r  series  of  the  form  (3)  yielding  roots  of  P(7(c),A)  with  k 
taking  the  values 

T]^'''ui^      (oT  j  —  1,  •  ■  ■ ,  r, 

for  some  tj  £  C  where  r;^/*"  is  the  principal  rth  root  of  rj  and  w  is  the  principal  rth  root  of 
unity.  An  immediate  consequence  of  this  type  of  sphUing  of  the  roots  of  P  is  that,  when 
evaluating  the  directional  derivative  (1)  for  the  spectral  abscissa  and  radius,  one  can  replace 
the  limit  infimum  in  the  definition  with  a  limit.  Moreover,  again  in  the  case  of  p  and  q,  this 
directional  derivative  always  exists,  possibly  taking  the  value  +00  but  never  the  value  —00. 
This  structure  also  indicates  why  the  spectral  radius  and  abscissa  are  Holder  continuous 
of  order  at  least  1/n.  Thus,  a  natural  extension  to  (1)  would  be  to  consider  Holder  type 
directional  derivatives  based  on  fractional  powers  of  e.  Although  it  is  possible  to  pursue 
such  an  approach  using  the  techniques  developed  in  Part  1  of  this  study,  we  do  not  do  so 
here. 


2     The  Directional  Derivative 

In  this  section,  we  apply  LemmEis  4  and  5  of  Part  I  of  this  study  to  obtain  information  about 
the  directional  derivative  defined  in  (1)  for  the  spectral  abscissa  and  radius.  We  continue 
to  assume  that  A:C  i— »  C"'*"  is  an  analytic  matrix  function  and  zq  £  C.  Also  let  A{zo) 
have  the  block  diagonal  decomposition 

Aizo)  =  A^°^  =  PDP-\     Z)=Diag(£)i,...,£>„),  (4) 

where  m  is  the  number  of  distinct  eigenvalues  of  A{zo)  and  the  kth  diagonal  block  Dt  is 
upper  triangular  with  constant  diagonal  At,  i.e. 

Dk  =  XkI+Nk, 

with  A'^it  strictly  upper  triangular  and  hence  nilpotent.  Let  n^  denote  the  order  of  Dj,  i.e., 
the  multiplicity  of  the  eigenvalue  At.  Suppose 


dz, 


{zo)  =  AY\  ioTl=l,--,s, 


and  let  Bn, . . . ,  Bmi  be  the  corresponding  diagonal  blocks  of  P    M(    P  for  /  =   l,--,s. 
With  this  notation  established,  we  have  the  following  results. 

Theorem  1   Choose  d  £  C .    If  for  some  k  £  A(zo)  ■=  {k  :  A^  €  S(zo)>   Re  At  =  0(20)} 
any  one  of  the  conditions 


Re  ^(tr  NkBt,)d,  <  0,  Im  ^^(tr  NkBk,)d,  =  0,  (5) 

/=i  (=1 

s 

Y^itTNi-'Bk,)di  =  0,j  =  S,...,nk,  (6) 


(=1 


15  violated,  then 

otherwise, 

where 


a''(-'o;rf)  =  +oo;  (7) 

a''(zo-d)>.^{zo,d),  (8) 


<p(zo,d)  —  maxi  ^^^^ :    k£A{zo)'>.  (9) 

Moreover,  if  conditions  (5)  and  (6)  hold  and  the  vectors 

{[iTNi-'Bki,--,  tT  Ni-'Bk,V  ■■   j=l,...,nk,  k£Ai(zo,d)} 

are  linearly  independent,   where  Ai(zo,d)  is  the  set  of  indices  m  A{zo)  which  achieve  the 
maximum  »n  (9),  then  equality  holds  in  (8). 


Partial  Proof  We  will  establish  the  first  statement  by  applying  Lemma  4  of  Part  I,  using 
the  correspondence  between  >4"'  in  Part  I  and  X]'=i  ^'^i  here.  It  follows  immediately 
that  (7)  holds  if  (5)  or  (6)  is  violated.  Indeed,  if  for  some  k  £  A(zo)  one  of  the  conditions 
(5)  or  (6)  fail  to  hold,  then  for  every  7  G  r{zo,d)  the  inequality 

Re(A(€)-  Xk)  <  6e  +  o(() 

fails  to  hold  for  every  (5  G  Tv,  for  some  branch  A(f )  of  the  splitting  of  At  along  the  curve  7. 
Consequently,  the  limit  defining  a'*(zo;d)  is  +00  for  every  7  G  r(zo,d)  and  so  a'^{zo,d)  = 
+00.  On  the  other  hand,  if  ^  >  Q'*{zo;d),  then  there  is  a  7  G  r{zo,d)  such  that 

a(y{())  -  a(zo) 
urn <  0, 

or  equivalently, 

Re(A(e)- At)  <  (5t  for  e  G  [0,(o], 

for  some  €0  >  0,  for  every  branch  A(e)  of  the  splitting  of  At  along  the  curve  7  ,and  for  every 
k  G  A(2o).  Now  by  applying  Lemma  4  of  Part  I  to  this  inequality  for  each  k  G  A{2o)  and 
then  letting  6  J.  Q'^{zQ\d),  we  obtain  inequality  (8).  The  proof  of  the  last  statement  requires 
the  implicit  function  theorem  and  is  given  in  [4].  n 

Theorein  2   Choose  d  £  C  and  set 

TZ{zo)  :=  {k  :  \\k\=  p{zo)]. 
We  will  consider  two  cases;  p(zo)  =  0  and  p{zo)  ^  0. 
1)  Let  us  assume  that  p{zo)  =  0  so  that  m  =  1  and  Ai  =  0.  //  any  one  of  the  conditions 

s 

J2(tT  N(-'Bu)d,  =  0,  j  =  2, . . . ,  Hi,  (10) 


(=1 


IS  violated,  then 

otherwise, 

where 


p''{zo;d)  =  +00  ; 
p''{zo;d)>azo,d),  (11) 


azo,d):=-\T{tTBi,)d,\ 
n    '-^ 


Moreover,  if  condition  (10)  holds  and  the  vectors 

{[trNi-'Bn,---,  UNr'BuV  :  j=l,...,n} 
are  linearly  independent,  then  equality  holds  m  (11). 


S)  Let  us  assume  (hat  p{zo)  ^  0.  If  for  some  k  G  1Z{zo)  any  one  of  the  conditions 

ReJt'^Yi^tT  Ni,Bk,)d,<0,  Im  A*' ^(tr  TVtBt,)'^/ =  0,  (12) 


^(trN^-'BH)(i/  =  0,  ;=3,...,rn,  (13) 


/=i 


is  violated,  then 
otherwise, 


p''{zo;d)  =  +00  ; 
p'{zo\d)>4izo,d),  (14) 


w, 


here 


»  t 


^(zo,d):=max{ --[Re  (A^  V (tr  Bi,)<i,)+  |V(tr  TVtBt,)^;!]    ■    k  en{zo)]. 

(15) 
Moreover,  if  conditions  (IS)  and  (13)  hold  and  the  vectors 

{[tr  Ni-'Btu--,  tTNl-'Bk,Y    :  j  =  l,...,nt,  keni{zo,d)} 

are  linearly  independent,  where  7vi(zo,d)  is  the  set  of  indices  in  7J(2o)  achieving  the 
maximum  in  (15),  then  equality  holds  in  (14)- 

Partial  Proof  First  consider  p(zo)  =  0.  It  follows  immediately  from  Lemma  5  of  Part  I 
that  (10)  must  hold  for  j  =  2  as  well  as  j  =  3,  •  •  -.n,  if  the  directional  derivative  is  to 
be  finite  valued.  Now  assume  that  (10)  holds  and  suppose  6  >  p^{zo;d).  Then  there  is  a 
7  G  T(zo,d)  such  that 

EA(0|<$;^|A(OI<n6e, 

where  the  sum  is  taken  over  all  branches  \{().  Letting  f)  [  0  yields  (11)  in  view  of  Lemma  3 
in  Part  L 

Next  assume  that  p(io)  ^  0.  The  proof  of  inequality  (14)  is  a  straight  forward  application 
of  Lemma  5  in  Part  I.  Observe  that  if  for  some  k  €  V-izo)  one  of  the  conditions  (12)  or  (13) 
fail  to  hold,  then  for  every  7  G  r(zo,d)  the  inequality 

\X{e)\  -  \Xk\<  6e  +  o{e) 

fails  to  hold  for  every  6  eTZ,  for  some  branch  X{e)  of  the  splitting  of  At  along  the  curve  7. 
Consequently,  the  limit  defining  p''{zo;d)  is  +00  for  every  7  €  T{zo,d)  and  so  p''{zo;d)  = 
+00.  On  the  other  hand,  if  6  >  p''(zo;d),  then  there  is  a  7  €  r(zo,d)  such  that 

,.     P(7(0)  -  P(zo)    .  c 
urn <  o, 


or  equivalently, 

|A(OI  -  |Ai|<6f  foree[0,eo], 

for  some  fo  >  0,  for  every  branch  A(f )  of  the  splitting  of  At  along  the  curve  7  ,and  for  every 
k  G  Tl{zo).  Now  by  applying  Lemma  5  of  Part  I  to  this  inequality  for  each  k  G  ^('o)  and 
then  letting  6  I  p'^{zo;d),  we  obtain  inequality  (14). 

The  proofs  of  equality  in  both  cases  makes  use  of  the  implicit  function  theorem  and  are 
given  in  [4].  q 

3     The  Nonderogatory  and  Semisimple  Cases 

In  this  section  we  consider  two  cases  that  illustrate  the  results  in  Theorems  1  and  2.  Let 
At  G  S(2o)  have  multiplicity  n^.  Then  A^  is  said  to  be  nonderogatory  if  it  possesses  only  a 
single  eigenvector.  The  nonderogatory  multiple  eigenvalues  are,  in  a  quantifiable  sense,  the 
most  generic  kind  of  multiple  eigenvalue  [1].  The  eigenvalue  At  is  said  to  be  semisimple  if 
it  is  associated  with  nj  independent  eigenvectors. 

Now  suppose  A{zo)  has  block  diagonal  form  as  described  in  the  previous  section.  Then 
At  is  semisimple  if  and  only  if  Nk  =  0,  and  At  is  nonderogatory  if  and  only  if  N^^~  -^  0. 
Therefore,  if  At  is  derogatory,  then  for  some  jo  <  "t  we  have  that  Nl  =  0  for  j  >  jo-  This 
observation  is  significant  in  the  interpretation  of  Theorems  1  and  2.  These  theorems  provide 
formulas  for  the  directional  derivatives  a  (zo\  d)  and  p''(zo;  d)  when  it  is  known  that  vectors 
given  by 

[tr;V^-'Bti,    ••,  tTNi-'Bk,V,  forj=l,...,nt 

are  linearly  independent.  However  these  vectors  cannot  be  linearly  independent  if  the  eigen- 
value At  is  derogatory  since  in  this  case  Nl~  =0  for  at  leeist  one  value  of  j  satisfying 
1  £  J  <  "i-  Thus  the  formulcis  that  are  provided  in  Theorems  1  and  2  are  at  best  sharp 
only  in  the  nonderogatory  case.  Such  an  outcome  could  have  been  predicted  from  the  out- 
set. As  has  been  indicated,  our  results  are  derived  from  results  for  polynomials  and  are 
extended  to  matrices  by  application  to  the  characteristic  polynomial.  Such  results  are  most 
applicable  to  nonderogatory  eigenvalues  since  the  multiplicity  of  these  eigenvalues  in  the 
characteristic  polynomial  agrees  with  their  multiplicity  in  the  minimal  polynomial.  In  this 
case,  the  characteristic  polynomial  accurately  reflects  the  behavior  of  the  eigenvalue.  This 
is  not  true  for  derogatory  eigenvalues. 

It  is  instructive  to  see  just  how  far  off  these  lower  bounds  can  be  by  considering  the 
semisimple  ceise.  Semisimple  and  nonderogatory  eigenvalues  lie  at  opposite  structural  ex- 
tremes whenever  n/c  >  1  (of  course,  if  nt  =  1,  then  the  eigenvalue  is  both  semisimple  and 
nonderogatory).  Thus,  we  would  expect  the  lower  bounds  for  the  directional  derivatives 
a''{zo;d)  and  p'^(zo,d)  to  be  most  inaccurate  in  this  case.  By  employing  the  resolvent 
theory  for  linear  operators  [5]  to  the  semisimple  case,  one  can  obtain  precise  formulas  for 
the  required  directional  derivatives.  This  is  done  in  [7,  Lemma  3.5].  As  an  immediate 
consequence  of  this  result  we  obtain  the  following  theorem. 

Theorem  3  Suppose  A(zo)  and  'K.(zo)  contain  only  semisimple  eigenvalues  of  A(zo).   Let 


d  £  C   and  denote  by  A'j.,,  for  I  =  l,...,nt,   the  eigenvalues  o/^'_,  Btidi  for  each  k  € 

JJ  The  spectral  abscissa:  The  ordinary  directional  derivative  Q'{zo,d)  exists  and  satisfies 
a'{zo;d)  =  a''{zQ-d)=     max       max  {Re  A^,}.  (16) 

2)  The  spectral  radius:  The  ordinary  directional  derivative  p'(zQ,d)  exists  and  satisfies 

{maxtg7e(^„)  maxi<K„,  \X'^.,\  if  p{zo)  =  0, 

-^  maxtgK(,„)  maxi<,<„^  Re  A^A;,       tf  p{zo)  7^  0. 

(17) 

Let  us  suppose  that  both  A{zo)  and  72(zo)  contain  the  single  element  k  and  that  the 
associated  eigenvalue  A^  is  semisimple.  Regarding  the  spectral  abscissa,  formula  (16)  says 
that  q'^{zo',  d)  is  the  maximum  of  the  real  parts  of  the  eigenvalues  of  the  matrix 

a 

Y^Btid,  (18) 

j=i 

whereas  the  lower  bound  for  a'^{zo;d)  given  by  the  right  hand  side  of  inequality  (8)  in 
Theorem  1  is  the  average  of  these  quantities.  This  behavior  is  repeated  in  the  case  of  the 
spectral  radius.  U  p(zo)  =  0,  then  p'^(zo;d)  is  the  spectral  radius  of  the  matrix  (18)  while 
the  lower  bound  given  by  (11)  is  the  modulus  of  the  average  of  the  eigenvalues  of  (18). 
Similarly,  if  p{zo)  ^  0,  then  p'^(zo',  d)  is  the  maximum  of  the  values 

ReXkX'u 


Pi^o) 


while  the  lower  bound  in  (14)  is  the  average  of  these  values  (since  Nk  —  0).  Apparently  the 
difference  between  the  actual  value  of  the  directional  derivative  (1)  and  the  lower  bounds 
given  by  Theorems  1  and  2  can  be  quite  great.  Nonetheless,  these  lower  bounds  indicate 
when  these  directional  derivatives  are  infinite-valued  and  in  the  generic  nonderogatory  case 
they  give  precise  values  for  the  directional  derivative. 

In  the  general  case,  when  the  eigenvalues  are  derogatory  but  not  semisimple,  little  precise 
information  is  known  beyond  the  general  lower  bounds  provided  by  Theorems  1  and  2.  The 
semisimple  case  illustrated  above  indicates  a  possible  upper  bound. 

4      The  Proximal  Normal  SubdifFerential 

In  this  section,  we  use  our  knowledge  of  the  directional  derivatives  a^{z;d)  and  p''{z;d) 
to  obtain  insight  into  the  structure  of  the  proximal  normal  subdifferential  for  the  spectral 


radius  and  spectral  abscissa.  A  vector  d  6  <C'  is  said  to  be  a  proximal  normal  subgradient 
for  the  function  f  :€'*->  TZ  At  zq  £  C  if 

,.      ■   ffjz)-  fi^o)-  {v,z-  Zq) 

liminf -„ TT >  -oo,  (19) 

where  the  inner  product  on  C  is  given  by 

(y,  w)  :=  Re  y^w. 

The  collection  of  all  proximal  normal  subgradients  at  zq  is  called  the  proximal  normal 
subdifferential  of  /  at  zq  and  is  denoted  by  d2f{zo).  Our  interest  in  understanding  the 
proximal  normal  subdifferential  for  the  functions  a  and  p  stems  from  the  rich  calculus 
that  is  available  for  this  subdifferential  and  its  limits.  For  more  information  about  this 
subdifferential,  we  refer  the  reader  to  Mordukhovich  [6]  and  Rockafellar  [8]. 

Let  Zq  €  C  and  suppose  that  A{zq)  has  the  block  diagonal  decomposition  described  in 
Section  2.  According  to  Theorem  1,  the  directional  derivative  a^(zQ\d)  is  bounded  below 
by  the  support  function  for  the  set 

1 


<I.(2o)     :=     {  —  [trSii, 


■,  ixB^,Y    ■■   keA(zo)  } 
+     [{(:[tr'NkBki,--,tTNkBk,V    ■■   ^€>l(2o),  ReC>0}U 

span{[trFr*Bu,-    •,  irNi'^Bk^V   ■  j  =  3,...,rn,  k  e  A(zo)  }], 
where  the  support  functional  for  a  set  Q  C  €'  is  defined  by  the  expression 


a^id)  :=sup(y,d). 
yen 

Similarly,  by  Theorem  2,  the  directional  derivative  p'*{zo;d)  is  bounded  below  by  the  support 
functional  for  the  set 

Eizo)     :=     {  —  [trS„,.--,  tr5i,r    :    ^G[0,2;r]} 
n 

+     span{[tr7v{''fln,--,  irN{~'Bi.Y   ■  i  =  2,...,n}, 
when  p(za)  =  0,  and  is  bounded  below  by  the  support  functional  for  the  set 

'    ;7d^[tr^*5H,--,trAr,B,,r 


*(^o)      :=      < 


+  ;77^[trS*i,---,  trBi,r 


nkp{2o) 


ken{zo,d),  ee[o,2w] 


+      [{(:^l[tTNkBku■■■^r^'kBk,]   :   k£ll{zo),ReC>0}U 


Tf7J-l 


span{[tr^i     Bti,--,triVi     Bk,Y    :  j  =  i, . .  .,nk,  k  €TZ{zo)  ]], 

when  p(zo)  ^  0.    The  vectors  defining  the  sets  $,  E,  and  *  are  our  prime  candidates  for 
proximal  normal  subgradients.  The  partial  results  that  we  have  been  able  to  obtain  to  date 


depend  upon  the  resolvent  theory  for  analytic  operator  valued  mappings  on  C  developed 
by  Baumgartel  [2].  The  key  result  that  we  require  in  this  regard  is  [2,  Theorem  2,  page  374]. 
In  order  to  understand  this  result,  one  must  first  understand  the  notion  of  a  A^-group  for 
Ai6S(^o). 

Given  Ajt  6  S(^o)  the  Ai-group  is  the  set  of  all  functions  X{z)  €  E(z)  such  that  X{z)  —> 
Ai  as  z  — ►  20-  The  A^-groups  are  well  defined  since  the  spectrum  E(z)  is  a  continuous 
multifunction.  The  projection  onto  the  invariant  subspace  for  the  entire  Ajt-group  is  called 
the  group  projection  for  A^.  The  following  remarkable  fact,  established  in  [2,  Theorem  2, 
page  374],  states  that  each  such  group  projection  is  locally  analytic. 

Theorem  4  Lei  zo  G  C.  Then  there  ts  a  neighborhood  Q  of  zq  »n  C  such  that  for  each 
Xk  €  l^l-o)  ''if  group  projection  n*:(-)  associated  with  A^  is  analytic  on  Q.  Moreover,  for 
each  Xk  G  5^(~o)  there  is  a  neighborhood  Gt  of  X^  m  (['  such  that 

Spec{A{z)nkiz))  =  ^z)nGi, 

for  all  z  £  ^  and  At  G  S(zo)  where  Spec(A{z)Uk{z))  denotes  the  spectrum  of  A{z)Uk{z). 

This  result  is  the  key  to  establishing  the  local  factorization  of  the  characteristic  poly- 
nomial used  in  Section  2  of  Part  I  of  this  study.  In  the  current  context,  we  simply  observe 
that  Theorem  4  implies  the  analyticity  of  the  trace  of  each  A^-group, 

niz)  :=  tr  A{z)nkiz),  (20) 

since  this  trace  is  simply  minus  the  coefficient  of  the  (A  —  A*)"*"^  term  in  the  characteristic 
polynomial  for  A(z)Pk{z).  From  Lemma  3  in  Part  I,  we  find  that  the  gradient  of  the  trace 
of  each  Afgroup  at  Zq  is  given  by 

Vrt(;o)  =  [trBu,--,  tr  S^,]^.  (21) 

By  combining  this  fact  with  the  observation  that  the  maximum  of  a  finite  set  of  real  numbers 
always  exceeds  (or  equals)  their  average,  we  are  able  to  establish  that  certain  members  of  the 
sets  ^(zq),  E(zo)  {'i(  p(zo)  =  0),  and  *(;o)  (if  p(-o)  ^  0)  are  proximal  normal  subgradients. 

Theorem  5  Lei  zq  G  C  and  suppose  thai  A(zo)  has  the  block  diagonal  decomposition 
described  in  Section  2.   Then 

—  [tr  Bu,  •  •  ■ ,  tr  'Bk>Y  €  d2a(zo),  (22) 

for  every  k  G  A{za). 

Proof  Let  k  G  >t(-o)-  Since  Tk{z)  is  analytic  near  ;o  and  a(z)  >  (l/nk)Tk{z)  for  z  near  zq 
and  k  £  Aizo),  we  have  that 

a(z)  -  a{zo)  -  —Re  Vnizofiz  -  zo)     >     —Re  [n(z)  -  n(zo)  -  Vn(zor(z  -  zo)] 


=     0(,\z-z,f). 


Thus,  by  (21) 


a(z)  -  0(20)  -  {^[U  Bku--,  tr  B^,],  z  -  Zo) 

liminf ;; ^ >  —00. 

[z-zoP 


D 
Theorem  6  Let  zq  €  C  and  suppose  that  A{zo)  has  the  block  diagonal  decomposition 
described  in  Section  2. 

1.  Suppose  that  p{zo)  =  0  so  that  E(io)  reduces  to  the  single  eigenvalue  A]  =  0.   Then 

^[trBn,--,  tTBuVed2p{zo),  (23) 

n 

for  every  ri  6  [0, 1]  and  0  6  [0,  27r]. 

2.  Suppose  that  p{zo)  ^  0.   Then 

— ^[tr  Sti,  •  •  ■ ,  tr  Bt.V  e  d2p(zo),  (24) 

for  each  k  G  7J(zo)- 

Proof    (1)  Let  rj  G  [0, 1]  and  9  G  [0,27r]  be  given.  Since  Ti{z)  is  analytic  near  zq,  we  have 

that 

77e'*       -  - 

p(z)-pizo)     -     {■^—[tTBtu--,iTBk,],z-zo) 
n 

>    -[\Mz)\  -  \n{zo)  +  Vn(zonz  -  zo)\] 

n 

=     Oilz-zof), 

which  establishes  the  result. 

(2)  Let  k  €  ■^(20)-  Since  Tt(z)  is  analytic  near  zq,  we  have  that 

/>(--) -P(^'o)     >     -[\tic{z)\  -  \nizo)\]  (25) 

\niz)\'  -  |rt(2o)P 

nk[\n(z)\  +  \n{zo)\] 

2Re  (rt(2o)(rt(2)  -  n{zo)))+  \n{z)  -  Tk(zo)\^ 

nk[\rk(z)\  +  \Tt(zo)\] 
2Re(n(2o)(rt(;)-r,(2o))) 

nk[\Tt(z)\  +  |n(zo)|] 
2Re  (rt(ro)(Vrt(2o)"(2  -  zp)  +  0(\\z  -  zpf))) 

M\rk{z)\  +  \Tt{zo)\] 
2Re  r,(zo)(Vnizonz-^  ^  ^^^^  _  ^^^^  ^^6) 


> 


2ni  |rit(2o)| 

'"\oilz-zon 


Re  XkVTt{zo)iz  -  zo)    ,  ^,,__      ,  ,,2 


nkP(zo} 

10 


where  the  inequality  (25)  follows  from  the  inequality  p{z)  >  (1/njt)  |rt(0)|  and  the  equation 
p(zo)  =  (l/nt)  ktC-o)!,  and  the  equality  (26)  follows  from  the  formula 

1  1  b-a 


(a  +  b)       26       26(a  +  6) 
and  the  fact  that  [\Tk{zo)\  -  \rk(^)\]  =  0(lz  -  2o|).  Therefore, 

p(z)-p(zo)     -     ( ^[tTBti,--,tTBk,].z-Zo) 

=     p(z)  -  p(zo) ^Re  (\k(^n{zoV(z  -  2o))) 

>     0(i~--''oP) 


Although  the  results  stated  in  Theorems  5  and  6  are  significant,  they  do  not  clarify  the 
status  of  the  remaining  vectors  in  the  sets  4>,  E  and  *.  We  have  yet  to  establish  that 
these  vectors  are  proximal  normal  subgradients.  Indeed,  we  are  uncertain  whether  or  not 
this  is  the  case.  The  difficulty  stems  from  the  fact  that  those  vectors  involving  the  matrix 
Nk  arise  from  the  bifurcation  of  the  eigenvalues  into  Puiseux-Newton  series  of  the  form  (3) 
where  the  power  p  satisfies  (1/n)  <  p  <  1-  Because  of  this,  the  individual  behavior  of 
these  eigenvalues  is  not  fully  captured  by  the  trace  of  their  associated  A^-group  since  these 
traces  are  analytic  at  zq.  However,  it  was  the  properties  of  these  traces,  especially  their 
analyticity,  that  allowed  us  to  establish  the  results  in  Theorems  5  and  6.  Consequently,  a 
different  approach  may  be  necessary  in  order  to  establish  whether  or  not  those  vectors  in 
the  sets  4>,  E  and  *  involving  the  matrix  Nk  are  proximal  normal  subgradients. 
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